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Using chiral anomaly, we discuss the pair creation of massless fermions under the effect of magnetic 
held B when an electric flux tube E parallel to B is switched on. The tube is axial symmetric and 
infinitely long. In the limit B > £, we can analytically obtain the spatial and temporal behaviors 
of the number density of the fermions, the azimuthal magnetic field generated by the fermions etc.. 
We hnd that the life time t c of the electric held is shorter as the width of the tube is narrower. 
Applying it to the glasma in high-energy heavy-ion collisions, we find that color electric field decays 
fast such as t c ~ QJ 1 with saturation momentum Q s . 
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PACS numbers: 12.38.-t, 24.85. +p, 12.38.Mh, 25.75.-q 
Schwinger mechanism, Chiral Anomaly, Color Glass Condensate 

Pair creation of charged particles in a classical strong electric field has been extensively discussedpj for long time 
since the discovery of Klain paradox, although such phenomena have not yet been observed. The particle creation is 



known as Schwinger mechanism^. Here, we are addressed with the problem in a quite different way from previous 
Qh| onesfH, 0, 0| and apply our results to color gauge fields (glasma) produced in high-energy heavy-ion collisions. 

The Schwinger mechanism have been discussed under spatially homegeneous or inhomogeneous ( and time dependent 
j" 1 1 1 or independent ) back ground electric fields with or without the effects of background magnetic fields. It has been 
treated in various formalisms, i.e. proper time method originally used by Schwinger, Heisenberg-Euler effective action 
formalism [31], or canonical formalism 4]. But only in a few papers back reactions of the produced particles have been 
discussed by taking account of the effects of background magnetic fields. In particular, there arc no such papers in 
which a magnetic field induced by the electric current of the particles is analyzed under the background magnetic 
field. On the other hand, the Schwinger mechanism in such a circumstance is quite relevant to the realistic cases such 
■^j- ■ as the glasma, neutron stars etc. Thus, it is an important issue to analyze the effects of the back reaction under the 
background magnetic fields. 

The glasma have recently received much attention. The gauge fields of the glasma are classical color electric 
and magnetic fields produced as initial states of gluons in high energy heavy ion collisions. We speculate that the 
decay of the glasma leads to thermalized quark gluon plasma (QGP) observed in recent RHIC experiments. Schwinger 
mechanism is a possible one for the decay of the color electric field. Indeed, such a mechanism was phenomenologically 
discussedQ about thirty years ago. Although the presence of the color electric field in the discussions was simply 
assumed, it has recently been confirmed on the basis of a fairly reliable effective theory of QCD at high energy, that 
is, a model of color glass condensate (CGC)@. 

According to a model of the CGC, the color electric and magnetic fields are homogeneous in the longitudinal 
direction and inhomogeneous in the transverse directions; we define that the longitudinal direction is parallel to the 
collision axis. Hence, we may think that they form color electric and magnetic flux tubes extending in the longitudinal 
direction. In our previous papers @, S Q, we have discussed a mechanism for the deca y o f the color magnetic field 
and have reproduced the instability of the glasma observed in numerical simulations [ill Irfl] . 

In this paper we discuss the decay of an electric flux tube under background magnetic field by taking account of 
the back reactions. The distinctive point in our paper is to use chiral anomaly without addressing the explicit forms 
of fermion's wavefunctions. Simply using the anomaly we can analytically obtain the number density of the fermions, 
the life time of the electric field, etc., all of which include the back reaction of the produced fermions in the electric 
field. Furthermore, we can analytically obtain an azimuthal magnetic field around the flux tube generated by the 
fermions. This is a remarkable feature in the use of chiral anomaly. However, the use of the chiral anomaly is limited 
only for the system with massless fermions under the presence of both electric and magnetic fields in the limit B ^S> E. 

Hereafter, we mainly discuss the pair creation of massless electrons and positrons in QED. The generalization to 
QCD would be straightforward. Before discussing the pair creation in an electric flux tube, we first explain how we 
use chiral anomaly in a simple model of a homogeneous electric field. We assume tentatively that the fields B and 
E are spatially homogeneous and parallel (or antiparallel) with each other. They are oriented in the z direction; 
B = (0, 0, B) with B > and E — (0, 0, E). Then, the energies of electrons with charge — e < and positrons with 
charge e > under the magnetic field B = \B\, but with E = are given by 
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E N = y/p 2 + 2NeB (parallel) and E N = + eB{2N + 1) (antiparallel), (1) 

where p z denotes momentum in the z direction and integer N > does Landau level. The term of "parallel" 
("antiparallel") implies magnetic moment parallel (antiparallel) to B. The magnetic moment of electrons (positrons) 
is antiparallel (parallel) to their spin. Thus, electrons (positrons) with spin antiparallel (parallel) to B can have 
zero energy states in the lowest Landau level; their energy spectrum is given by -Ejv=o = \Pz\ > 0. They are 
"massless" states. On the other hand, the other states cannot be zero energy states; their energy spectra are given 
by \Jp1 + eBM > eBM > eB with positive integer M. They are "massive" states whose masses increase with B. 

When the electric field is switched on, the pair creation of electrons and positrons with the energies En begin to 
occur. However, it is not probable that any states are produced with an equal production rate. Indeed, high energy 
states are hard to be produced while low energy states are easy to be produced. Thus, the "massive" states are hard 
to be produced by weak electric field \E\ -C B. In particular, they cannot be produced in the limit of B — > oo. 
Only products in the limit are the pairs of electrons and positrons in the "massless" states. We assume such a limit 
B ^S> \E\ in our discussion below. 

Note that owing to the Fermi statistics only fermions with p z = are produced, since the other states with p z ^ 
have already been occupied. After the production of the fermions with p z = 0, their momenta increase with time, 
p z (t) — ±e J dt'E, owing to the acceleration by the electric field. Hence we obtain the momentum distribution fl2| 
n(p z ) oc 9(pF(t) — p z )9(p z ) for positrons with charge e > and n(p z ) oc 9(pp(t) + p z )9(—p z ) for electrons. The Fermi 
momentum pp is given by PF(t) — e J Q dt'E(t'). Here we have assumed that the electric field is parallel to B and is 
switched on at t = 0. 

We note that electrons move to the direction antiparallel to E while positrons move to the direction parallel to E. 
Therefore, both electrons and positrons created by the electric field have right handed helicity when E being parallel 
B, while they have left handed helicity when E being antiparallel B. The helicity of the particles never change in the 
course of their propagation. 

The key point in our discussion is to use the equation of chiral anomaly, 

d t (n R -n L ) = ^E(t)B (2) 

where np), ( ul ) denotes number density of right ( left ) handed chiral fermions; nn,L = ("I , 7o(l±75) , I')/2 in which the 
expectation value is taken by using a state of electrons and positrons created in the electric field. Here, we have used 
spatial homogeneity of the chiral current j'5, that is, diyjs = 0. Since we have assumed sufficiently strong magnetic 
field B E, the back reaction of the produced particles on B is negligible. ( Although the magnetic moments of the 
fermions partially screen the back ground field B, the effect is negligible in the limit of B 3> E. ) 

In accordance with the chiral anomaly, the rate of chirality change is given by the product of E and B. It is 
important to note that the chirality is identical to the helicity in the case of massless fermions. Thus, the anomaly 
equation describes the rate of helicity change. As the helicity of the particles does not change in the course of 
their propagation, the equation describes the rate of the particle production or annihiration. Since only particles 
with right ( left ) handed helicity are produced in the limit B ^> E when E is parallel to B (E being antiparallel 
to B ), riR = 2n and nj, = with the number density n of electrons when E parallel to B, while tir — and 
np = 2n when E antiparallel to B. Thus, we find that the pair production is determined by the anomaly equation 
2d t n(t) = e 2 E(t)B/47r 2 with E{t) > 0. 

Similarly, the chiral anomaly describes the pair annihilation dtn = e 2 E (t) B / Air 2 < when E(t) < 0. The annihi- 
lation may occur subsequent to the pair creation. Indeed, after the particle production the electric field E > loses 
its energy and vanishes because the field accelerates the charged particles. Then, the field changes its sign, i.e. E < 
since a Maxwell equation dtE(t) = — J with non vanishing electric current J = 2en(t) of the particles, dictates that E 
should not remain zero ( see later. ) Negative E in the anomaly equation implies the decrease of the number density 
n with time, which is caused by pair annihilation. 

In order to discuss back reactions of the particles, we need to consider an energy conservation between the energies 
e(t) of the particles and the energy of the field, 

dt (e{t)+ l -E 2 {t))=Q (3) 

where we have neglected the contribution of magnetic field induced by the electric current of electrons and positrons. 
This is allowed in the system with spatially homogeneous fields E and B. Later we take into account the contribution 
when we discuss an electric flux tube. 
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The energy density e(t) of the particles are given such as e(t) = 2n(t)pp(t) /2 = n(t)pp(t). This originates from the 
fact that the momentum distribution h(p z ) ( n(t) = ± J ±o ° dp z h(p z ) ) is given by h(p z ) = n 9(pp(t) — \p x \)9(\p z \) as 
we have discussed. For example, the energy density of electrons is given such that J_ dp z h(p z )\p z \ — npp/2. 

Now we have three equations to solve for obtaining the pair production rate dtn(t), etc., 

e 2 / 1 \ 

2d t n(t) = —E(t)B, dt(e(t) + ^E 2 (t)j = 0, and e(t) = n(t)p F (t) (4) 

with pp(t) = J Q dt'eE(t'). It is easy to solve the eqQ with initial conditions E(t = 0) = Eq > and n(t = 0) = by 
assuming magnetic field B independent on t. 



E(t)=E Q casU—t) and n(t) = =L=^ (5) 

TT 2 f^E 



with a = e 2 /4ir. 

In Fig. 1 we have shown the behaviors E(t) — E cos(irt/2t c ) and n(t) = aE °2 tc \ sm(nt/2t c )\ with Eo = 1, t c = 1 
and a ~ 1/137, where t c = The formula holds rigorously in the limit of B Eq where only right handed 

or left handed particles are produced. The oscillating features of E(t) and n(t) have agreed with those obtained in 
the previous papers [l|, in which the wavefunctions of electrons under the effects of the electric and magnetic fields 
have been explicitly used. We also note that the screening of the magnetic field B by the magnetic moments of the 
fermions is negligible in the limit B 3> E since the number density n of the fermions is proportional to \/~B. 

Here we wish to explain why the number density n and the electric field E oscillate. When the electric field 
E > is switched on at t — 0, the pair production of the particles with right handed helicity begins to occur and 
their number density increases. Owing to the acceleration of the particles by the electric field, the Fermi momentum 
pp(t) — L dt'eE(t') and their energies increases. On the other hand, the energy of the electric field gradually decreases 
and vanishes at t = t c . At the time the number density takes the maximum value. Then, the electric field changes 
its direction ( E < for t > t c ) and becomes strong with time, as discussed before. Accordingly, the number density 
decreases after t = t c since 2d t n(t) — -£^jE(t)B < 0. This decrease is caused by the pair annihilation of electrons and 
positrons. Since the direction of the electric field changes after t c , the direction of the particle acceleration also change 
so that the Fermi momentum p F (t) begins to decrease. Consequently, electrons and positrons are moved to overlap 
with each other so that the pair annihilation may occur to make the number density decrease. Eventually, all of the 
particles with right handed helicity vanishes at t = 2t c . At the time the strength of the field E(t = 2t c ) = —Eq < 
becomes maximum. Since E is negative and n = Q at t = 2t c , a new pair creation of the particles with left handed 
helicity begins to occur after t > 2t c . Their number density increases in accordance with the chiral anomaly, 
2dtn(t) — —j^E(t)B > 0. Therefore, the oscillation of E(t) and positiveness of n{t) as shown in Fig.l arises. 

Here we wish to make a comment why the electric field increases after t = t c . Since all the positrons ( electrons ) 
move to the direction of B ( —B ) at t = t c , the electric field anti-parallel to their electric current becomes strong by 
gaining its energy from the particles, once such an electric field is present even if it is small, on the other hand, the 
electric field parallel to the current vanishes losing its energy. Therefore, the electric field increases with its direction 
anti-parallel to B when t > t c . 

In the above discussion, we have assumed that the mass m e of electrons and positrons is vanishingly small. The 
assumption may hold as far as E > m 2 .. Thus, when the electric field becomes small such as E < m 2 by losing the 
energy, the anomaly in eq((2]) cannot be used. Actually the pair creation can hardly arise in such a case. However, 
the field loses its energy due to the acceleration of the particles. Thus, it vanishes ( E = ) and increases with the 
opposite direction ( E < ) in acoordance with the Maxwell equation as explained above. Once it arises with its 
direction opposite to that of the electric current, it gains the energy. Consequently, the field can becomes strong such 
as E > m 2 so that the assumption of the masslessness approximately holds again. Although the detail in the real 
behavior of E is different with that in Fig.l when E is small, the feature of the oscillation we found still remains. 

It is important to note that although the formula of the chiral anomaly involves all of the quantum effects, the 
mechanism why the number density increases or decreases is not explicit in the formula. However, we understand 
that the increase or decrease of the number density originates with the pair creation or annihilation. This is because 
the results shown in this paper have agreed with those obtained by calculations]!] explicitly including the quantum 
effects of the pair creation or annihilation. 



Now we discuss the pair production in an axial symmetric electric flux tube, namely the electric field extending 



4 



FIG. 1: electric field E(t) (solid line) and number density 
n(t) (dashed line) 

infinitely in the z direction, but has finite width R in the transverse directions. In this case the quantities of E, n 
and pf depend on the cylindrical coordinate r — \/ x 2 + y 2 and time coordinate t. ( Since the wavefunctions in the 
Lowest Landau level are given by Z m exp(— eB\Z\ 2 /A) with Z = x + iy and integer m > 0, the charged particles 
can be located within the flux tube for sufficiently large B such as R \j\feB. ) Thus, the formula of the chiral 
anomaly takes the same form as the one in eqQ, although the electric field E(r,t) and the number density n(r,t) 
depend on r and t. ( In the equation of the chiral anomaly, we may take divja = d z j z + ^ i=x y dijf = 0, because 
of the homogeneity in the z direction and the fact that J2i=xy^df = Yli=x y ^{^757^) — f° r the "massless" 
states in the Lowest Landau level. ) Similarly, the energy density e takes the form e(r,t) = n(r,t)pF(r,t) as eq(j4|. 
in which the Fermi momentum is given by pf(j~, t) = J Q dt'eE(r,t'). We also take into account Maxwell equations 
such as d t Be{r, t) — d r E(r, t) and d t E{r, t) = MljEillifl) _ j( r) where J(r, t) represents the current of electrons and 
positrons flowing in the z direction. Bg(r,t) denotes azimuthal magnetic field generated by the current. Therefore, 
we have the following equations to solve, 

e 2 

2d t n(r, t) = -^E{r, t)B, e(t) = n(t)p F (t), (6) 



and Maxwell equations, 



d t Bg(r,t) = d r E{r,t), d t E{r, t) = d ^ rBe ^ _ J( r ,t). (7) 



where B is assumed independent of r and t because back reactions of the fermions on B are small in the limit B ^> E. 

In order to solve the equations we must represent the current J(r, t) in terms of n(r, t) and pf{t, t). It can be done 
by imposing the condition of energy conservation, 

dt J d 3 x ^ 2 M) + £ 9 2 M) +e(7 . t) ) = j d 3 x(E(r,t)dtE(r,t) + Bg(r,t)d t B e (r,t) + d t e{r,t)) =0. (8) 

The condition implies that the energy of the electric field is transmitted into the energies of the particles and the 
azimuthal magnetic field Bg. 

Using the equations in eq© and eq([7]), we rewrite the condition in eq© such that 



J d 3 x(E(8 t E - drfr-Bg) ) + g te ^ = J d 3 x (_ JE + = J d 3 xE(en + ^Bp F - J) = 0, 



(9) 



where we have performed the partial integration in r. The energy conservation must hold for any initial conditions of 
E(r, t = 0) = E (r). Therefore, we obtain 

e 2 

J[r, t) = en{r, t) + ^Pf{t, t)B = 2en(r, t), (10) 
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FIG. 2: electric field £(r) at t = 0.2 (small dots) 
and E(r) at t — 0.6 (large dots) 
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FIG. 3: number density n(r) at t — 0.2 (small dots) 
and n(r) at t — 0.6 (large dots) 



since ^^pp(r,t)B — #^ J". dt'eE{r,t')B — en(r,t) with n(r, < = 0) = 0. It states that both of electrons with charge 
density —en and positrons with en constitute the electric current J since the velocity of electrons is —1, while that of 
positrons is +1 in the unit of light velocity c = 1. This formula of J — 2en has been used in the discussion just above 
eq©. 

Using J(r, t) in eq (fT0]) as well as equations in eq© and eq©, we derive the equation of motion of the electric field, 



d 2 E(r,t) = (d 2 



1 e s B 



(11) 



where we have used the formula dfJ(r,t) — -£^BE(r,t). 

We solve the equation (TrTj) with the initial conditions n(r,t = 0) = 0, Bg(r,t — 0) = and E(r,t — 0) = 
Eq exp(— r 2 /R 2 ). The initial conditions implies that the electric flux tube with radius R is switched on at t = 
when any particles and magnetic field Bg are absent. Using the initial conditions E(r, t = 0) = Eq exp(— r 2 / R 2 ) and 
dtE(r, t = 0) = derived from Maxwell equations (J7J), we first obtain the solution E in eq (fTTj) . and then we derive n 
and Bg, 



E(r,t) 
n(r, t) 
B e (r,t) 



kdk cos(t^k 2 + m 2 )J (kr) exp(-k 2 R 2 /4), 



e 2 BE R 2 



,, sm(tx/ k 2 + m 2 ) ,. , , , o „ 9 , , 
fcdfc — v , \ . ^o(H exp(-fc 2 i? 2 / 4 ) 



y/k 2 + m 2 



16tt 2 

E oR 2 r, 2 „ sin(tVA: 2 + m 2 ) 



k 2 dk- 



\Jk 2 + m? 



Ji(kr)exp{-k 2 R 2 /4), 



(12) 
(13) 
(14) 



with m 2 = -£^B, where Ji(kr) denotes Bessel functions. Since the integration over k is dominated by small k in 
the limit R — ► oo, we find that the solutions in ea fT^)) ~ ea (rT4l are reduced to the previous solutions in eq© in the 
limit. We note that the effect of Bg on the particles is negligible compared with the background magnetic field B 
since Bg oc 1/ \J~B in the limit B — > oo. 

We have shown the temporal and spatial behaviors of electric field E(r, t) in Fig. 2, number density n(r, t) in Fig. 3, 



and azimuthal magnetic field Bg in Fig. 4, respectively, where the parameters m = 1, EqR /2 = 1, e 



= 1 and 



i? = 1 have been used; the scale r = 10 in the figures corresponds to the width R of the electric flux tube. Small 
and large dots represent the behaviors at t = 0.2 and t = 0.6, respectively. We find that the stronger electric field 
produces the particles more. It leads to the fact that the number density of the particles is larger as r is smaller. 
Thus, the electric field E(r, t) vanishes faster as r is smaller since E(r, t = 0) is stronger as r is smaller. We also 
find that the life time of the electric field with finite R is shorter than that of the field with R — oo. ( The life time 
may be defined roughly as E(r — R,t — t c ) — 0, although the field oscillates with time. ) This is because typical 
momentum k dominating in the above integration is on the order of R . Thus, the life time is approximately given 
by t c ~ (Vm 2 + R- 2 )~ 1 . 

The quantities E(r,t), n(r,t), and Bg(r,t) oscillate with time owing to the similar reason explained above. In 
particular, electric field E{r, t) never disappear by completely losing its energy. This originates with the fact that the 
produced particles are almost free and do not interact with externals e.g. external heat bath so that their energies 
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FIG. 4: azimuthal magnetic field B${r) at t = 0.2 
(small dots) and Bg(r) at t — 0.6 (large dots) 

never dissipate. On the other hand, when we take into account their interactions with external heat bath and use 
fermionic momentum distributions at finite temperature, the electric field would decay and disappear because the 
current J would vanish owing to the dissipation of the particle's momenta in the heat bath. 

Finally we wish to make a comment on the decay of the glasma produced in high-energy heavy- ion collisions. The 
color electric and magnetic fields generated in the early stage of the collisions form flux tubes and their field strenghts 
are almost identical to each other. Although we have assumed that the gauge fields are Abelian and magnetic field 
B is not a tube-like but homogeneous and B 3> E, we may approximately apply our results to the analysis of the 
decay of the color electric field. In particular we use the result that the life time t c of electric field is approximately 
given by (\/R~ 2 + m 2 ) -1 . In the glasma, m 2 is of order a s gB where gB ~ Q 2 S with saturation momentum Q s ( 
a s = g 2 /47r ~ 1/47T with gauge coupling constant g ~ 1 ), while the width of the flux tube R is given by Qj 1 . Hence, 
we find that the color electric field decays rapidly such as t c ~ Qj 1 ~ 0.1 fm/c ~ 0.2 fm/c for Q s = 1 GeV ~ 2 GeV. 
The smallness of the width is the main origin of the rapid decay. Although our result does not include expanding effects 
of the glasma, the result suggests that the color electric field decays sufficiently fast as required phenomenologicallv [T5| . 
In addition, the result is consistent with our assumption of the independence of B on time. Indeed, B decays very 
slowly [Tol [lij and does not vary so much at least within the period of Q^ 1 . 

To summarize, using the chiral anomaly we have discussed the pair creation of massless charged fermions in an axial 
symmetric electric flux tube under the effect of strong magnetic field. We have analytically obtained the spatial and 
temporal behaviors of their number density, the electric field and the azimuthal magnetic field induced by the current 
of the fermions. Our calculations have been performed without addressing explicit forms of particle's wavefunctions. 
These quantities of the free massless fermions are exact in the limit of B 3> E. Applying our results to the glasma 
in high energy heavy- ion collisions, we have found that the color electric field decay sufficiently fast to be consistent 
with QGP phenomenology. 

The author express thanks to Drs. H. Fujii of University of Tokyo and Dr. K. Itakura of KEK for their useful 
discussion and comments. He also thanks to Dr. Tanji of University of Tokyo for his numerical result of the momentum 
distribution when m e — > 0. 
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